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Abstract 

We propose a new method for constructing the quasi-exactly solv- 
able (QES) potentials with two known eigenstates using supersymmet- 
ric quantum mechanics. General expression for QES potentials with 
explicitly known energy levels and wave functions of ground state and 
excited state are obtained. Examples of new QES potentials are con- 
sidered. 
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1 Introduction 

A potential is said to be quasi-exactly solvable (QES) if a finite number 
of eigenstates of corresponding Schrodinger operator can be found exactly 
in explicit form. The first examples of QES potentials were given in [1-4]. 
Subsequently several methods for generating QES potentials were worked out 
and as a result many QES potentials were found [5-13] (see also review book 
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[14]). Three different methods based respectively on a polynomial ansatz 
for wave functions, point canonical transformation, supersymmetric (SUSY) 
quantum mechanics are described in [12]. 

Recently, an anti-isospectral transformation called also as duality trans- 
formation was introduced in [15] . This transformation relates the energy 
levels and wave functions of two QES potentials. In [16] a new QES poten- 
tial was discovered using this anti-isospectral transformation . 

The SUSY method for construction of QES potentials was used in [10- 
12]. The starting point of this method is some initial QES potential with 
n + 1 known eigenstates. Then applying the technique of SUSY quantum 
mechanics (for review of SUSY quantum mechanics see [17, 18]) one can cal- 
culate the supersymmetric partner of the QES potential. From the properties 
of the unbroken SUSY it follows that the supersymmetric partner is a new 
QES potential with n known eigenstates. 

In addition SUSY was used to develop some generalized method for the 
construction of the so-called conditionally exactly solvable (CES) potentials 
in [19, 20]. The CES potential is the one for which the eigenvalues problem 
for the corresponding Hamiltonian is exactly solvable only when the poten- 
tial parameters obey certain conditions. Such a class of potentials was first 
considered in [21]. 

In the present paper we develop a new SUSY method for generating QES 
potentials which in contrast to the previous papers [10-12] does not require 
the knowledge of the initial QES potential for the generation of a new QES 
one. 

2 Supersymmetric quantum mechanics 

In the Witten's model of supersymmetric quantum mechanics the SUSY 
partner Hamiltonians H± read 




+ V ± (x), 



(1) 



where 
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V ± {x) = l -(w\x)±W>{x)), W '(x) = ^^-, (3) 

W(x) is referred to as a superpotential. In this paper we shall consider the 
systems on the full real line — oo < x < oo. 

Consider the equation for the energy spectrum 

H ± ^(x) = E±^(x), n = 0,1,2,.... (4) 

As a consequence of SUSY the Hamiltonians H + and H- have the same 
energy spectrum except for the zero energy ground state. The latter exists 
in the case of the unbroken SUSY. Only one of the Hamiltonians H± has 
a square integrable eigenfunction corresponding to zero energy. We shall 
use the convention that the zero energy eigenstate belongs to H^. Due to 
the factorization of the Hamiltonians H± (see ([!])) the ground state for H 
satisfies the equation B~ip$ (x) = the solution of which is 

^(x) = C ~ exp (- J W(x)dx ] j , (5) 

Cq is the normalization constant. Here and below C denotes the normal- 
ization constant of the corresponding wave function. From the condition 
of square integrability of wave function ipQ (x) it follows that superpotential 
must satisfy the condition 

sign(W(±oo)) = ±1. (6) 

Note that this is the condition of the existence of unbroken SUSY. 

The eigenvalues and eigenfunctions of the Hamiltonians H + and if_ are 
related by SUSY transformations 

E n+1 = E t > E = °> ( 7 ) 

VWi(z) = -±=B + i;+(x), 1>+(x) = - 7 L=B-^ 1 {x). (8) 

V E n V E n+1 

For a detailed description of SUSY quantum mechanics and its application 
for the exact calculation of eigenstates of Hamiltonians see reviews [[17], |1S[[ . 
The properties of the unbroken SUSY quantum mechanics reflected in SUSY 
transformation (|7]), (H) are used for exact calculation of the energy spectrum 
and wave functions. In the present paper we use these properties for the 
generation of the QES potentials with the two known eigenstates. 
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3 QES potentials with the two known eigen- 
states 



Suppose we study a Hamiltonian H_, whose ground state is given by (|5|). 
Let us consider the SUSY partner of H_, i.e. the Hamiltonian H + . If we 
calculate the ground state of H + we immediately find the first excited state 
of H- using the SUSY transformation (|7|), (H). In order to calculate the 
ground state of H + let us rewrite it in the following form 

H + = H W + e = B+Bf + e, e > 0, (9) 

which leads to the following relation between potential energies 

V+(x) = V W (x) + e, (10) 

where e is the energy of the ground state of H + since we suppose that 
has zero energy ground state, Bf and V_ (x) are given by (|) and (D with 
the superpotential W\(x). 

As we see from (||) the ground state wave function of H + is also the ground 
state wave function of H_ and it satisfies the equation B^ip^x) = 0. The 
solution of this equation is 

^+(x) = C+ exp (- J W 1 {x)dx S j , (11) 

where for square integrability of this function the superpotential W\ (x) satis- 
fies the same condition as W(x) @. Using (0) and (§) we obtain the energy 
level Ei = e and the wave function of the first excited state ip\{x) for 
From (]To|) we obtain the following relation between W(x) and W\(x) 

W 2 (x) + W\x) = W?(x) - W[{x) + 2e. (12) 

Previously the same equation was used in the case of the so-called shape 
invariant potentials to obtain the exact solutions of Schrodinger equation 
(see also reviews [0, 0). We consider a more general case and do not 
restrict ourselves to the shape invariant potentials. Note, that (O) is the 
Riccati equation which can not generally be solved exactly with respect to 
W(x) for a given Wi(x) and vice versa. 
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The basic idea of this paper consists of finding such a pair of W(x) and 
Wi(x) that satisfies equation dl2|) . It has been recently suggested by us in 
P5fl . For this purpose let us rewrite equation (JfJ) in the following form 

W' + {x) = W-(x)W+(x) + 2e, (13) 

where 

W + (x) = W 1 (x) + W(x), (14) 
W_(x) = W!(x)-W(x). 

This new equation (|13|) can be easily solved with respect to W-(x) for a 
given arbitrary function W + (x) or with respect to W + (x) for a given arbitrary 
function W-(x). Then from (H) we obtain superpotentials and Wi(rr) 



which satisfy equation (12). 



3.1 Solution with respect to W-(x) 

In this subsection we construct the QES potentials using the solution of 



equation (|1^) with respect to W-(x) 

W-(x) = (W'+(x) - 2e)/W+(x), (15) 

where W + (x) is some function of x. Note, that the superpotentials W(x) 
and W\(x) must satisfy condition @. Then as one may see from (|T4] ) W+(x) 
must satisfy the same condition @ as VF(a;) and do. 

Let us consider continuous functions W + (x). Because W + (x) satisfies 
condition (|J) the function W + (x) must pass through zeros. Then as we see 
from flT5D iy_(x), and thus W(x), W\(x) have poles. In order to construct 
the superpotential free of singularities suppose that W+(x) has only one 
zero at x = x$ with the following behaviour in the vicinity of Xq W + (x) = 
WAxq){x — Xq). In this case the pole of W-{x) at x = x can be cancelled 
by choosing 

e = W' + {x )/2. (16) 
Then the superpotentials free of singularities are 

W{x) = I (w + (x) - (W'(x) - W'(x ))/W + (x)) , (17) 
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Wi(x) = \ (w+(x) + (W[(x) - Wl(x ))/W+ 



, x 
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Substituting the obtained result for W(x) into (^) we obtain QES poten- 
tial V- (x) with explicitly known wave function of ground state (||) and wave 
function of the first excited state. The latter can be calculated using (|TT|) 
and © 

$-{x) = C- W+(x) exp (- J Wi(x)dx) . (18) 

It is indeed the wave function of first excited state because W+(x) has one 
zero. 

We may choose various functions W + (x) and obtain as a result various 
QES potentials. The functions W+(x) must be such that ipo(x) and ipx(x) 
are square integrable. If the eigenfunctions ipo(x) and ipx(x) belong to the 
Hilbert space of square integrable functions in which the Hamiltonian is 
Hermitian then these functions must be orthogonal 



<V l^i >=-C C 1 



exp ^— J dxW + (x] 



0. (19) 



The wave functions must also satisfy appropriate boundary conditions. 

To conclude this subsection let us consider explicit examples. Choosing 
W + (x) = A (sinh(o;x) — sinh(axo)) we obtain the well known QES potential 
derived in ||] by the method elaborated in the quantum theory of spin 
systems. Note that the case xo = corresponds to Razavy potential [||]. 
Following the proposed SUSY method this example is considered in details 
in our earlier paper [^5| . 

Consider the function W + (x) in the polynomial form 

W + (x) = ax + bx 3 , a>0, 6>0 (20) 
which gives a new QES potential 

T , . , 1 . o . q ab & b 2 (■ 3ab 3b a , . 

VJx) = -(a 2 - 126 s 2 + —x 4 + —x 6 + — + — — - - -. 21 

v ; 8 V ; 4 8 8(a + bx 2 ) 2 8(a + bx 2 ) 4 v ; 

The energy levels of ground and first excited states are Eq = 0, E± = a/2. 
Note, that two energy levels of this potential do not depend on the parameter 
b. The wave functions of those states read 

^j-( x ) = C (a + &x 2 )3 /4 e- x2(2a+bx,2)/8 , (22) 

$-{x) = Cfx(a + bx 2 ) 1/4 e- x2(2a+bx2)/8 . (23) 

It is worth to stress that the case 6 = correspond to linear harmonic 
oscillator. 
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3.2 Solution with respect to W + (x) 

Equation (|13|) is the first order differential equation with respect to W+(x). 
A general solution can be written in the following form 



W+(x) = exp (^J dxW-(x)j 2e J dxexp (- J dxW-(x)j + A 



(24) 



here A is the constant of integration. 

In order to simplify solution ( p4f) let us choose W-(x) to be of the form 



W„(x) = -<f>"(x)/<f/(x), 
and suppose that 4>'{x) > 0. Then 

W+(x) = (2e0(x) + A)/0'(x). 



(25) 



(26) 



Note that the constant A can be included into the function <p(x) and thus 
for W+(x) we obtain 

W+(x) = 2e0(x)/0'(x). (27) 
Finally for superpotentials W(x) and W±(x) we have 



^(x) = (^"(x) + e0(x))/0'(x), 
W 1 {x) = {^<P'\x)-e<P(x))/<P\x). 



(28) 
(29) 



Using this result for wave functions of the ground state with the energy 
Eq = and excited state with E± = e we obtain 

^-{x) = C o -(0'(x))" 1/2 exp \-e J dx<f){x)/(j)'{xf) , Eq = 0, (30) 
^-(x) = C^(f)(x)((f)'(x)y 1/2 exp f-e J dx<j>{x) / <j>' {x)\ , E^ = e, 

where function 4>(x) must be such that these wave functions are square inte- 
grable. The condition of orthogonality in this case can by written similarly 
to© 



exp 



dx<f>(x) / <f>' (x) 



0. 



(31) 
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QES potential VL(x) is given by ([|) with superpotential (|28|) . Choosing 
different <p(x) we obtain different QES potentials with explicitly known two 
eigenstates. We shall consider a nonsingular monotonic function <f>{x) with 
one node. Then ipi(x) also has one node and thus corresponds to the first 
excited state. 

In conclusion of this subsection let us consider an explicit example. Let 
us put 

<f>{x) = ax + bx 3 /3, a, 6 > 0. (32) 

Note, that the case 6 = corresponds to a linear harmonic oscillator. The 
function ( |3~2"D generates the following superpotentials 

W(x) = ex/3 + (b + 2ae/3) X , - , (33) 

a + bx A 

W x (x) = ex/3 + (-6 + 2ae/3) ^ - (34) 

a + bx z 

which as we see satisfy condition @. 

Substituting W(x) into (|3|) we obtain the following QES potential V-(x) 
and its SUSY partner V+(x) 



v -^-^ x2 + ^ + J^y +R " (35) 



where 



A_ = A + = e 2 /9, B_ = b + lae, = R+ = -^-(36 + 4ae), 

3 186 

D_ = — ^-(27a6 2 + 24a 2 6e + 4a 3 e 2 ), D + = -\(9ab 2 - 4a 3 e 2 ). 

186 v ; ' + 186 v ; 



Using ( p0| ) we obtain the wave functions of the ground and first excited states 

^j-( x ) = Cq (a + bx 2 y 1/2 - at/3b exp(-ex 2 /6), (37) 
^-( x ) = C- (ax + bx 3 /3){a + bx 2 )- 1/2 ~ ae/3b exp(-ex 2 /6). (38) 



The same result for QES potential fl35|) and wave functions 
can be obtained using the method described in section 3.1 and taking the 
function W + (x) to be of the form (p7|). 
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It is interesting to stress that in the special case e = 3b/ 2a the QES 
potential V-(x) reads 

VLfs) = -^x 2 + — - — — + — (39) 

v ; 8a 2 a + bx 2 (a + bx 2 ) 2 4a v ; 

and can be solved exactly. To see this note that in this special case the super- 
potential Wi(x) = ex/3 corresponds to superpotential of a linear harmonic 
oscillator. Then (x) and, as a result of (|10]), V+(x) are the potential 
energies of the linear harmonic oscillator 

The fact that V+(x) corresponds to the linear harmonic oscillator follows also 
directly from fl36D because the coefficient D + in the considered case is equal 
to zero. Therefore in this case H + is the Hamiltonian of the linear harmonic 
oscillator and we know all its eigenfunctions in explicit form. Using SUSY 
transformations (0), (Q) we can easily calculate the energy levels and the 
wave functions of all the excited states of H-. Note that in this special case 
V-(x) can be treated as CES potential and it corresponds to the one studied 



in na 1201, 23, 24 



As far as we know the potential in general form ( p5|) has not been pre- 
viously discussed in the literature. This potential is interesting from that 
point of view that in the case of e = 3b /2a this potential is the CES one for 
which the whole energy spectrum and the corresponding eigenfunctions can 
be calculated in the explicit form. 
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